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Tremendous progress has been made towards the solution of the binary-black-hole problem in 
numerical relativity. The waveforms produced by numerical relativity will play a role in gravitational 
wave detection as either test-beds for analytic template banks or as template banks themselves. As 
the parameter space explored by numerical relativity expands, the importance of quantifying the 
effect that each parameter has on first the detection of gravitational waves and then the parameter 
estimation of their sources increases. In light of this, we present a study of equal-mass, spinning 
binary-black-hole evolutions through matched filtering techniques commonly used in data analysis. 
We study how the match between two numerical waveforms varies with numerical resolution, initial 
angular momentum of the black holes and the inclination angle between the source and the detector. 
This study is limited by the fact that the spinning black- hole- binaries are oriented axially and 
the waveforms only contain approximately two and a half orbits before merger. We find that for 
detection purposes, spinning black holes require the inclusion of the higher harmonics in addition to 
the dominant mode, a condition that becomes more important as the black-hole-spins increase. In 
addition, we conduct a preliminary investigation of how well a template of fixed spin and inclination 
angle can detect target templates of arbitrary spin and inclination for the axial case considered here. 



I. INTRODUCTION 

Initial LIGO has reached its design sensitivity [H 
and, along with other ground-based detectors around the 
world, is taking science data. The matched filtering tech- 
nique employed to search for gravitational wave signals 
from inspiraling and merging compact binaries is at its 
best when accurate representations of the signal are used. 
One of the best sources for early detection is the inspiral 
and merger of two black holes. The accurate waveforms 
of the binary black hole (BBH) signal come from two 
sources: post-Newtonian theory and numerical relativity. 
Waveforms from post-Newtonian theory are well under- 
stood during the inspiral phase of the binary evolution 
up to an as yet unknown point close to the final plunge of 
the black holes. With numerical relativity, we attempt to 
determine the point at which the approximation breaks 
down and to supply the waveform for the last orbit (s) 
and merger. Because black holes in the 30 — SOOM© Q 
mass range are expected to merge at frequencies that are 
in the most sensitive part of the LIGO frequency band, 
the computation of BBH waveforms has been a prior- 
ity for the numerical relativity community. Fortunately, 
remarkable progress continues to be made in numerical 
relativity, especially in solving the last orbits and merger 
of binary black hole systems. 

Since the breakthroughs in numerical relativity that 
allowed the evolution of BBHs through an orbit, merger 
and ringdown [3,, ^, ^ , many groups have achieved similar 
success. The numerical relativity has now begun evolv- 
ing many orbits 0, 0j B Q and explorin g th e parameter 
space such as mass ratio and spins [lOl I 111 12, 13, 14, 
m, m, [13, [H, [H, Each evolution of a BBH system 
requires large amounts of computational resources. The 



resource requirements for longer evolutions with unequal 
masses and spins are larger still. Numerical relativity ef- 
forts will benefit from a priori knowledge of where best 
to spend those resources to create an appropriate library 
of waveforms. Additionally, the cost of doing data anal- 
ysis for a large template bank with many parameters is 
also computationally expensive. Knowledge of the sen- 
sitivity of the waveforms to the parameters will guide 
efforts in data analysis and numerical relativity. As the 
library of numerical relativity waveforms grows, we can 
determine how the waveforms will be best employed in 
the search for gravitational waves and the characteri- 
zation of their sources. They may be used in conjunc- 
tion with post-Newtonian waveforms through stitching, 
as templates for binary mergers, or to test the quality of 
template spaces built from post-Newtonian and/or ap- 
proximate waveforms. 



In [2l| , a comparison of the closeness between current 
gravitational wave search template banks and numerical 
relativity waveforms from non-spinning BBHs has been 
performed. The numerical relativity waveforms were 
used as target signals against which the template banks 
currently in use by LIGO to search for inspiral gravita- 
tional waves were tested. They found matches greater 
than 0.96 for many of the analytic and approximate tem- 
plate families at 10— 120A/q. In 9], a phenomenological 
family of waveforms was proposed to model the coales- 
cence of the BBHs using a hybrid method that combines 
analytical and numerical relativity waveforms. They 
achieve matches greater than 0.99 at 3OA'/0 — 130Mq for 
non-spinning BBHs. 

In this paper, we conduct a first-step investigation into 
the ramifications on matched filtering of including spin 
in the BBH waveforms provided by numerical relativ- 



2 



ity. The set of spinning BBH waveforms that we use are 
provided by the Penn State group and were pubhshed 
in [iBl in addition to several non-spinning BBH wave- 
forms. The spinning binaries are axial configurations 
with the initial black holes having spins of equal magni- 
tude where one is aligned and one anti-aligned with the 
orbital angular momentum for a set of four spin param- 
eters, a = J/A'P = {0.2,0.4,0.6,0.8}. The physical pa- 
rameters of the waveforms are the numerical resolution, 
A, and the initial angular momentum of the black holes 
in terms of the spin parameter, a. Because the configu- 
ration is axial, the intrinsic parameters of the waveforms 
do not include the angles between the spins and the axis 
of the orbital angular momentum. We also include the 
radiation modes given by i and m of the spin weighted 
s = —2, spherical harmonics, -2Yim, which are the in- 
clination and azimuthal angles between the source and 
detector in the source frame. These waveforms contain 
approximately two and a half orbits. Because of this, we 
restrict the total mass range to to > 50Mq, set by the 
initial frequency of the evolution and our choice of how 
we enter the LIGO band. 

Both numerical accuracy, e.g. truncation errors, and 
astrophysical accuracy, e.g. initial data choices, will play 
a role in determining the viability of numerical relativ- 
ity BBH waveforms acting as templates. Requirements 
may be more stringent when characterizing the sources 
of gravitational waves; however, in this paper we focus 
on the use of numerical relativity waveforms as potential 
templates for detection, not for parameter estimation. 
Requirements for detection were first placed on wave- 
forms from BBHs generated by numerical relativity in 
references 0, i2^] . This early work preceded the success- 
ful solution of the BBH problem by many years, but acts 
as a guide for determining the constraint on numerical 
resolution accurate enough for data analysis purposes. A 
similar method was employed more recently in reference 
[2^ in connection with BBH evolutions of equal-mass, 
non-spinning black holes over several orbits. Ref. [23j . 
made a prediction of a maximum match that numerical 
waveforms will resolve. In addition, we did a preliminary 
study of the impact that numerical errors can have on the 
faithfulness of numerical waveforms using a Zerilli-based 
toy model in [24] . We now employ similar tests as [1| and 
[231 and verify the predicted behavior of the match with 
resolution for our BBH evolutions. 

Often in numerical relativity, waveforms are extracted 
in terms of coefficients of _2^m the dominant mode be- 
ing the quadrupole mode {£ = 2, |to| = 2). For com- 
pact binary inspiral searches, restricted post-Newtonian 
templates are commonly used for detection [25]. These 
templates include only the dominant harmonic in the am- 
plitude while including as much information as possible 
in the phase since phasing is the more important issue in 
matched filtering. Corrected-amplitude templates have 
been considered in [2^ and found to reduce the signal- 
to-noise ratios for LIGO and add features to the detection 
and parameter estimation for Advanced LIGO [13, . 



In [23] , it was found that using the £ = 2 mode was was 
good enough for detecting gravitational waves of non- 
spinning, equal mass binaries. We include in our analysis 
the angle from the spherical harmonics, 9 and 0, where 
9 is related to the inclination angle of the binary with 
respect to the detector but in the source frame. We find 
that, in general, we will need to include higher modes to 
accurately represent spinning BBH waveforms. We also 
explore the sensitivity of the match with the spin param- 
eter of the black-holes' initial angular momentum and 
how well a reduced template bank would do in matching 
with a target template of arbitrary spin and inclination. 

The paper is organized as follows. In fjll] we explain 
the techniques we have employed in our BBH code and 
in our data analysis algorithms used in the rest of the 
paper. We present our results in S jllll including a de- 
scription of the equal-mass, spinning and non-spinning 
evolutions that produced the waveforms under consider- 
ation including their spectra. We first assess the quality 
of our waveforms in terms of resolution through an anal- 
ysis of the convergence properties of our BBH code in 
mil Al We also present how the match statistic can re- 
flect the convergence properties of our waveform in ijlllBI 
In mil CI we present our study of the spinning templates 
in terms of varying inclination angle. Finally, in mil Dl 
we look at using a specific numerical relativity waveform 
as the template and check its ability to match a target 
of unknown spin and inclination. We conclude in ijlVl 

II. TECHNIQUES 

A. Numerical Relativity Techniques 

We solve the fully nonlinear Einstein equations in the 
3-1-1 form. Bowen-York |2^] puncture [30] initial data 
representing the two black holes is specified on a spacelike 
slice t = 0. This allows us to specify the initial mass, 
linear momentum, coordinate location and spin of each 
black hole. The BSSN [Ml, [H, jp formulation of the 
Einstein equations in the x form [5| is used to evolve this 
data forward in time. The gauge conditions are modified 
versions of the l-|-log lapse and F-driver shifts, 

doa = -2aK, (1) 
doP' = 3/AB\ (2) 
doB' = dor-vB\ (3) 

where do = dt—P^dj. We choose = 2 as in [3^. The ad- 
vection term [3^ djV^ removes certain zero-speed modes of 
the system as analyzed in [s^. These gauge choices, with 
the exception of this advection term, were found to be 
important for long-term stable and accurate evolutions 
of head-on collisions without excision (36j . The punc- 
ture locations are determined by integrating the equation 

dtx\t) = (}\x{t)) [13. 

The Weyl scalar ^'4 = Rap-ysn'^'nT'^n^'fn^ measures out- 
going gravitational radiation at large distances from the 
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source. The orthonormal tetrad {n", l", m", m"} is con- 
structed from the spherical coordinate basis vectors as in 
[ssl ]. We decompose \l/4 into _2y«m on spheres of constant 
coordinate radius. 

We use the Cactus infrastructure for paralleHzation, 
I/O and parameter handhng. The initial data is com- 
puted using the TwoPunctures [39] thorn, a spectral code 
which solves the momentum constraint for the conformal 
factor of our conformally flat spatial metric. We use the 
moving punctures approach without excision 4,^ so the 
singularity is not treated specially, apart from ensuring 
that it does not lie on a grid point initially. The evolu- 
tion is performed using our BSSN thorn which is auto- 
matically generated using the Kranc code generation 
package. ^'4 is also computed using a Kranc-generated 
thorn. For mesh refinement we use Carpet |4l|, and our 
Movelt infrastructure for specifying the grid structure 
based on the locations of the black holes. 

The code uses finite differences, and all spatial deriva- 
tives are computed using centered fourth order stencils, 
apart from the shift advection terms which are computed 
using upwind, lop-sided, fourth order stencils. This gives 
increased accuracy over using centered stencils through- 
out. The evolution stencil thus has three points adjacent 
to any given point in all three directions. For time inte- 
gration we use fourth order Runge-Kutta with a Courant 
factor dt/dx' = 1/2. 

We use 9 levels of box-in-box mesh refinement, where 
the outermost (base) grid covers the domain G 
[—320,320]. This allows the outer boundary to be 
causally disconnected from the computation of '^4 at 
r = 30m until well after the wave has passed. The next 
three grids cover the domains a;' G [—160, 160], [—80,80], 
and [—40,40]. These grids are all fixed in place through- 
out the simulation. The remaining 5 levels generically 
contain two grids each centered on one of the black holes. 
These grids are cubical of half-side 10, 5, 2.5, 1.25 and 
0.625. Our moving grids infrastructure handles the case 
of the two grids on a given refinement level overlapping 
by replacing them with a single grid which is the smallest 
grid containing both of them. Each grid has half the grid 
spacing of its parent grid, and we typically refer to the 
resolution of the finest grid as ft./. This grid structure 
allows the extraction sphere dX r — 30m to not intersect 
any refinement boundaries. 

We use Berger-Oliger mesh refinement as implemented 
by Carpet, and we refine in time as well as in space. 
The evolution on the refined grids requires a bound- 
ary condition provided by the coarse grid. This is im- 
plemented by adding additional points outside the re- 
finement boundary, which are filled by interpolation 
from the coarser grid. Since the evolution stencil is 
of size 3, this boundary must consist of at least three 
points. Since the fine grid must be updated twice as 
frequently as the coarser grid, the interpolation is in 
time as well as space. As noted in [4l|, for evolution 
equations which are first order in time but second order 
in space, this scheme leads to instabilities. One solu- 



tion is to interpolate from the coarse grid only at the 
start of the Runge-Kutta sub-stepping algorithm, not 
at every sub-step. Since on each sub-step the boundary 
points are not computed, additional buffer points are re- 
quired outside the finer grid. These points are also filled 
by interpolation at the start of the sub-stepping algo- 
rithm. The total number of boundary points required is 
stenciLsize X time-integrator _sub-steps = 12. We use fifth 
order interpolation in space, and second order interpola- 
tion in time. This lowers the overall accuracy expected 
from 4th to 2nd, but using higher order time interpola- 
tion is prohibitively expensive in terms of computational 
resources. We find that using 9 boundary points instead 
of 12 does not affect the results significantly, but it re- 
duces the computational cost [13]. We expect that the 
second order error from the mesh refinement boundaries 
is small, and that at low resolutions the scheme will ap- 
pear fourth order accurate. 

The initial data parameters are given in Tab. HI The 
non-spinning Rl parameters are taken from [43. b ased 
on the quasi-circular sequence found in |43|. |44 E^. 
The spinning runs are those reported in [l5|. These 
runs consist of two equal-mass black holes whose ini- 
tial spins are oriented such that one is aligned with 
the orbital angular momentum and the other is anti- 
aligned, and the magnitudes of the spins are equal. The 
runs are labeled as S0.05-S0.20 and have spin parameters 
a — {0.2, 0.4, 0.6, 0.8}. The black holes are located at po- 
sitions (0,±?/, 0), have linear momentum (=pP, 0, 0), spin 
(0, 0, ±5*), and bare puncture masses m±. The irreduce- 
able masses (measured from the areas of the apparent 
horizons) are toi = m2 = to/2 where to = toi + m2 is 
the total mass of the binary. 

TABLE I: Initial data parameters 



Model y/m P/m a m+/m m-/m mAum/m Jadm/tti^ 

Rl 3.257 0.133 0.483 0.483 0.996 0.868 

S0.05 2.95 0.13983 0.2 0.4683 0.4685 0.98445 0.825 

SO. 10 2.98 0.13842 0.4 0.4436 0.4438 0.98455 0.825 

S0.15 3.05 0.13547 0.6 0.3951 0.3953 0.98473 0.825 

S0.20 3.15 0.13095 0.8 0.2968 0.2970 0.98499 0.825 



B. Data Analysis Techniques 

Our numerical waveforms are calculated during evolu- 
tion in terms of the Newman-Penrose ^'4, where 

rm^'4(t,r) = ^-2Q„(i,r)_2>lm(^,'^), (4) 

and r is the extraction radius. The relationship between 
^1^4 and h^{t) and hx{t) is given by 

^A{t)^^{h+{t)-^h^,{t)). (5) 
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To carry out the matched filtering analysis, we compute 
h+{f) and h^{f) directly from the Fourier transform of 
the real part of ^'4 such that 



(6) 



This avoids issues regarding integration constants that 
arise from the time domain integration of ^4 Our 
initial data contains gravitational radiation which does 
not correspond to that present in an astrophysical situ- 
ation. This radiates away and is visible in the waveform 
as an initial data pulse. We remove this pulse from each 
waveform. 

Our analysis of the match of the numerical relativity 
waveforms will follow along the lines of the matched fil- 
tering procedure for detecting gravitational waves. That 
is, given two time domain waveforms hi{t) and h2{t), 
the scalar or inner product between these two functions 
is defined as 



(/ii|/i2) =4Re 



frr. 



hiif)Kif) 
Shif) 



(7) 



where the domain [/min, /max] is determined by the detec- 
tor bandwidth and h{f) stands for the Fourier transform 
of the respective time series. Sh{f) denotes the noise 
spectrum for which we use the initial LIGO noise curve. 
The fact that our numerical waveforms only contain a 
few orbits before merger means that they do not span 
the entire LIGO frequency band. The initial orbital fre- 
quency of these runs varies depending on the value of the 
spin between approximately 0.016/m and 0.024/m. The 
most stringent lower limit on /min would be 0.024/m. 
We impose the condition that the signal-to-noise ratio 
would be coming entirely from the domain spanned by 
our numerical waveforms such that /min is calculated for 
each mass of the template. This fixes the lowest mass for 
which the templates could be useful to SOM©. 

The match statistic [13], is properly defined as the 
maximized overlap between the signal and the template, 
however, we will use it as a measure between two tem- 
plates hi and /i2. There are two extrinsic parameters: 
the time of arrival of the signal, to, and the initial phase 
of the orbit when it enters the LIGO band, $. First, con- 
sider the maximized overlap between templates defined 
with only maximization over time given by 



lax [hi, /l2] 
0[hi,h2] 



maxO[ft-i, ft-2] , 

to 

(fel|fe2) 
V{hl[hl){h2[h2 



(8) 
(9) 



In Eq. ([8]), the maximization is to be understood as the 
maximum overlap between hi and /12 obtained by shifting 
the template in time, i.e. h{t) — > h{t -\- to) leading to the 
numerator being written in frequency domain as 



max (/ii 1/12} 

to 



max 4Re 

to 



^i(/)/i*(/)e2"/*o 

/min 



Skif) 



df. 

(10) 



The Fourier transform is replaced by a discrete fast 
Fourier transform, transforming the integral into a dis- 
crete sum. When computing the match without phase 
optimization, we will typically compute the match of the 
templates using only 

We also compute the match with an optimization over 
the phase, $ of the template in addition to to- A wave- 
form of arbitrary initial phase $ is written as 



/»(/) = /»+(/) cos $ + /ix(/) sin $. 



(11) 



Given an arbitrary waveform, hi, the phase optimization 
over the template, /i2, can be carried out using the nor- 



malized templates 



'i+j X I 



where i runs 



over 1,2. In reference to the match we will always be 
referring to frequency domain templates and for ease of 
notation we shall refer to h as h in the match formulae. 
The typical match fi^ is given by 



typ = max maxO[/ii+, ft.2] 

to '^2 



~ max VO[ei+, e2+]2 + 0[ei+, e2x]^ , (12) 

to 

where we have assumed that 62+ and 62 x are nearly or- 
thogonal, i.e. (e2+|e2x) ~ 0. This approximation is valid 
only for angles near {9 = 0, (p = 0) where the contri- 
bution from the modes other than ^ = |m| = 2 is zero 
because of the spherical harmonic values. For a general 
waveform at nonzero 9 and (j), the overlap of the corre- 
sponding h+ and hx is of the order of a few percent. In 
our spinning waveforms, we find a maximum deviation 
from orthogonality of ^ 3% at an angle that approaches 
tt/2. In order to avoid uncertainties in the matches at 
large angles, we construct an orthonormal waveform ba- 
sis as done by [48i, |49i] and outlined below. 

Any arbitrary polarization other than x) can be 
expressed as a linear combination of these linearly in- 
dependent basis vectors. The phase optimization of the 
match can then be done in terms of the new orthonormal 
basis vectors, one of which we choose to be ei+ and the 
other ei± given as 

ei± = (eix - ei+(e^+|eix>) (1 - {ei+\eixfy^^^ . (13) 

One can see that (ei+jeij,) = by construction. Given 
such orthonormalized basis vectors for two sets of pa- 
rameters i = 1,2, we can calculate the upper and lower 
bounds on the phase optimized match [50|. 

In this new orthonormal basis, we rewrite the typi- 
cal match from Eq. (|12p by matching the -|- polarization 
of one template while optimizing over the phase of the 
other. This mimics the situation in which one of the 
waveforms acts as the template and can be maximized 
over its phase, $2, while keeping the phase of the second 
template, $1, fixed, as follows 



typ 



max maxO[/ii+, ^,2] 



to * 



= max VO[ei+, e2+]2 -I- 0[ei+,e2±]^ ■ (14) 

to 
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In general, the phase of the signal could take any value, 
and one would like to know the best and worst possible 
values of the match. The expression for the best match 
(the upper bound) is given by [HOl 



Mbest = max max maxO[/ii, /i2] 

to *i *2 



(15) 



max 

to 



A + B 



A-B 



in which the phases of each template are optimized. The 
minimax match is given by the case when one maximizes 
the phase of one of the templates but minimizes over the 
other. This is to mimic the worst case scenario when 
the signal phase is such that it gives lower matches even 
when maximized over the template phase, given by [soj ] 



M, 



minimax 



maxminmaxO[/ii, /12] 

to ^2 *1 



(16) 



max 

to 



A + B 



A~B 



and for both cases the functionals A, B, C are written 
in terms of the orthonormal basis functions (ei+,eij_) 
and (e2+, 62^) corresponding to the two templates being 
compared 

A = (ei+|e2+)^ + (ei+|e2±)^ , 
B = (eij^|e2+)^ + (ei^|e2±)^ , 
C = (ei+|e2+) (ei_L|e2+} + (ei+|e2_L) (ei_L|e2±) . 

Note the symmetry of the formulae upon interchange 
of the two templates which arises from the fact that only 
the relative phase between the two templates should mat- 
ter. Because of this, which template is maximized and 
which minimized is interchangeable. We will refer to the 
target's phase being minimized and the template's phase 
being maximized since this makes sense in a detection 
scenario. 



III. RESULTS 

We present our results in three categories: the varia- 
tion of the match with resolution, including convergence 
tests of the waveforms; the variation of the match with 
inclination and spin; and how well we can differentiate 
between a target of unknown inclination angle and spin 
in the mass range of 50Mq to 3OOM0 given a template 
of fixed parameters. The spectra of the BBH waveforms 
from the models Rl and S0.05-S0.20 are shown in Fig. [2 
in which we plot \h{f)\. The 6 = lines correspond to 
the £ = \m\ =2 mode. The final black hole in all these 
runs settled to the same final spin with a spin parameter 
a « 0.66. 

In Fig. [H we can see that when 9 = 0, the spectra for 
each value of a become very similar for / > 0.05m. This 




f[m] 



FIG. 1: We plot the spectrum of the wave versus frequency 
of waveforms from different initial spins in a log- log plot. 
The spinning cases all have resolutions of m/40 and the non- 
spinning, m/38.4. We present three cases of inclination angle, 
9, for each set of spins, and the units are /code = /phys/m. 



is because all of the configurations we evolved settled 
down to very similar final black holes. When 9 deviates 
from zero, however, we can already see from the spectra 
that we get more variation between the different spinning 
waveforms. This is in part due to the larger resolution 
requirements for the higher harmonics, however, as we 
shall see in the following results, the resolution differences 
do not account for all of the variation. 



Firstly, in ^III Al we investigate the convergence of the 
dominant, £ = \m\ = 2, mode for non-spinning BBHs, 
the convergence of the spinning waveforms was published 
in [1^ to be between third and fourth order. In mil Bl 
we investigate the quality of the numerical waveforms as 
templates in matched filtering in terms of the resolution 
including spin and inclination angle. Next, in mil CI we 
show how the inclusion of modes £ < 5 affects the match 
calculations when compared to just using the dominant 
mode waveform £ = |m| = 2 as a function of the spin. 
Additionally, we determine the effect of using a finite 
extraction radius on the quality of the matches. Lastly, 
in mil Dl we compare the different spin waveforms to each 
other, including the dependence on 9 but keep the masses 
the same. This is equivalent to checking the faithfulness 
[50| of the waveforms assuming the orientation of the 
source for the different sources to be the same. While 
this assumption will need to be relaxed in a fully general 
treatment for data analysis, it suffices to demonstrate the 
importance and impact of using higher than dominant 
modes in the analysis of these merger waveforms. 
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A. Numerical Convergence 

A convergence test of the spinning series of runs 
called S0.05-S0.20 is given in Here we investi- 

gate of the quality of our numerical waveforms with 
a convergence test of the equal-mass, non-spinning 
Rl series of runs. These runs are names Rla-Rlf 
and only differ in their resolution. The suffix corre- 
sponds to the grid spacing of the finest grid surround- 
ing each black hole. The grid spacings are Aa-/ = 
{to/25, to/32, to/38.4, to/44.8, to/51.2, to/57.6}. We 
consider the convergence properties of the £ = to = 2 
mode of ^4 computed on the coordinate sphere at r = 30. 
In Fig. [3] we plot Re[^'4^], and Fig. [2] shows that ^-4^ 
computed at r = 30 appears to converge monotonically 
with resolution to a continuum solution, and that the 
high resolution results agree well with each other. Simi- 
lar results (not shown) are obtained for arg ^'4'^ and the 
coordinate locations of the black holes. Using the runs 
Rla, Rib, and Rlf we demonstrate fourth order conver- 
gence of the amplitude and phase of ^"4'^ in Fig. d] and 
Fig. [5l Were we to plot only the highest resolution runs, 
the fourth order convergence would be lost; i.e. there is 
a source of error which spoils the convergence at high 
resolutions. However, as can be seen from Fig. [2l the 
waveforms still appear to converge to a continuum solu- 
tion, just at a different order. 



0.0026 


1 


1 


1 


1 1 










Rla 


0.0025 








Rib - 


0.0024 








-Rlc 










\ Rid 


0.0023 








,\\ Ri«? - 


0.0022 


// 
// 






\\\Rif\ 

V 


0.0021 


/' '' 
// / 






\ \ 
\\ ''. \ 

V '■ * 


0.002 


/ 1; / 









195 200 205 210 215 220 
t 





003 




0.002 - 




0.001 - 












0.000 - 






Pi 






-0.001 - 




-0.002 - 




-0.003 - 




50 



100 



150 
t 



200 250 300 



olutions. The feature at f = 40 is due to the gravitational 
radiation present in the initial data, and will be cut out of 
the waveform before taking Fourier transforms for data anal- 
ysis. 
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FIG. 4: Convergence of |'I'4' | using the highest and the two 
lowest resolutions. C3, C4, C5 are the scaling factors expected 
for third, fourth and fifth order convergence respectively. We 
see that a — h matches most closely with C4 x{b — f) indicating 
fourth order convergence. 



FIG. 2: Monotonic change in |^'4'^| as resolution is increased. 
It appears that the function is converging to a continuum 
solution. 



B. Variation of Templates with Resolution 

We can follow the numerical convergence properties 
through matched filtering since we know how the numer- 
ical codes approximate the exact solution to the partial 
differential equations discretely. For finite differencing 



this is expressed as 

h{t)^ho{t) + c{t)AP + 0{AP+') (17) 

where ho{t) is the exact solution found when the resolu- 
tion goes to infinity, A is the grid spacing (A = Ax = 
Ay = Az oc At), p is the order of the truncation error 
and c(t) is a time dependent scaling independent of the 
grid spacing. Since our code is approximately fourth or- 
der accurate, we expect to scale with p — A. One can 
calculate hoit) using Richardson extrapolation of three 
discrete, convergent solutions of different resolution. We 
have computed ho(t) for the Rl series of runs using Rla, 
Rib, and Rlf and found that ho{t) coincides to within 
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FIG. 5: Convergence of arg*I'4'^. As in Fig. lU we see fourth 
order convergence. 

10~^ in amplitude with the waveform evolved with the 
finest resolution, A/ = m/57.6. In light of this, we will 
use our finest resolution run, indicated as /i(A/), in place 
of the Fourier transform of ho{t). 

We can predict the behavior of the maximized overlap 
as a function of resolution between two waveforms that 
differ only in their resolution by expanding the match 
equation, Eq.®, about A^* — 0. Note that when A^' is a 
constant in time, we can express Ea. (|17p . in the Fourier 
space as hi = hg + cA^, where i is running over each 
waveform. We expand the match for the case of two 
templates at two different resolutions to be as general as 
possible, i.e. 0[/ii(Ai), /?,2(Aj)] as follows 

+ 0(3)(Af). (18) 

The above equation was expanded for two matches of 
different resolution but the match with the Richardson 
extrapolated solution can be recovered by setting A^ = 
up to the order expressed. This equation indicates that 
the mismatch, (1 — Omax), goes like A^^. Fig. ([6]) shows 
the dependence of the Omax on the resolution for three 
choices of the total mass, m = {50, 100, 2OO}Ajf0. The 
least squares fit is done for several possible values of A'^p 
including p = 3, 4, 5. The best fit is found for p = 4, and 
is the fit pictured in Fig.®. 

Because of the computational cost of generating solu- 
tions to the Einstein equations at three resolutions for ev- 
ery waveform of interest, Flanagan and Hughes [2^ sug- 
gested computing the match between templates of neigh- 
boring resolutions, such as Omax[^(Aa), h( A.i,) ] for exam- 
ple. A variation, used by Baumgarte et al [23|, is to com- 
pute the series of matches as Omax[/i(Ai), /i(A/)], where 
the first template runs over all the resolutions available, 
denoted by A^, and the second template's resolution is 
fixed to the finest available. 
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FIG. 6: The maximized overlap, Omax [/ii, ft/], from Eq. (|8]) 
is plotted as a function of grid size for three total masses 
indicating fourth order convergence. The data are given by 
the points and the least square fitting by the lines. The p = 4 
fit differs for different masses because of the weighting by the 
current LIGO noise curve Sh{f)- 



We plot the Flanagan and Hughes type of overlaps, 
Omax[/i+(Aa),/i+(Ab)] and Omax[/i-H(Af,), /i+(Ac)] and 
three cases where we vary the resolution of one template 
while keeping the second template fixed to the finest 
resolution, Omax[/i+(Ai), /i+(A/)] where i = a,b,c, in 
Fig. ([7]). These matches are unoptimized over the phase, 
but optimized over the time to and follow the definition 
given in Eq. ([5]). 

The threshold of match above which the resolutions 
will be "good enough" for detection is typically set to 
0.98. Clearly from the figure we can conclude for the = 
0, a resolution of A;, = m/32.0 was enough for the entire 
mass range to reach a minimum match of 0.98 and that, 
for a lower resolution of A^ — m/25.6, we only achieved 
Omax > 0.98 for total masses larger than IOOMq. One 
of the main features of the matches with resolution is 
the decrease in the match toward smaller masses. This is 
partly due to the low number of gravitational wave cycles 
we have present in the waveforms. The numerical errors, 
however, will increase as we increase the total number 
of cycles so the choice of resolution will depend on the 
number of orbits evolved. As we evolve more cycles, we 
will also be extending the total mass range into smaller 
masses. 

We now have to choose which match to use in assessing 
the quality of the templates. The matches are expected 
to get closer to one when phase optimization is included. 
Since the phase of a given signal waveform depends on the 
detector orientation and several such variables, it could 
differ from the template up to a constant phase factor and 
the match would vary depending on what the phase is. 
The worst case would occur when the arbitrary phase is 
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FIG. 7: The maximized overlaps, Omax[/ii+, /i2+], are plotted 
as a function of mass for several resolutions of our equal- 
mass BBH series of runs. Labels (a, b, c) stand for the coarse 
runs Rla, b, c with resolutions m/25.6, m/32.0, m/38.4, while 
Rlf has resolution m/57.6. The overlap between adjacent 
resolutions Omax['i+(Aa), h+{At) has deceptively high values 
compared to the matches with the finest resolution. 



such that there is maximum destructive interference be- 
tween the target and the template, and sets up the lower 
bound on the possible matches, A^minimax- If one does 
not add any phase factor to the numerically generated 
template , one gets what is called the typical match and 
if the phase factor is chosen for maximum constructive 
interference, it yields the "best" match, Afbost- The value 
of Mbcst is a measure of the closest distance between two 
templates. We'll use the minimax match for the rest of 
the paper to report our results in terms of lower bounds 
and use the typical and best matches to illustrate the 
range over the different type of matches. 

To demonstrate the differences in the match with 
choice of optimization, we tabulate examples of opti- 
mization in Tab, im for few selected total masses between 
the coarsest and finest runs. No optimization at all is 
listed in the second column, optimization over just Iq 
from Eq. ([5]) in the third, the typical match Eq. (|14p in 
the fourth, minimax, Eq. (|17p . in the fifth and finally 
best match, Eq. I|16p . in the final column. The phase 
optimized matches, the last three columns, exceed the 
the threshold of Omax > 0.98 for even the least resolved 
waveform when matched with the finest resolution wave- 
form. 

We also calculated the minimax matches, see Eq. ([T7|. 
for all the Rl resolutions, although we do not plot them 

> 



minimax 



here, and found that for all the resolutions Af, 
0.99 in the mass range under consideration and show the 
same variation as the overlap. These results hold only 



TABLE IL Comparing matches with different optimizations 
between the coarsest (Rla) and the finest (Rlf) runs for se- 
lected total masses. The second column is the normalized 
overlap between the two waveforms without optimization over 
to, the third column is the match with optimization over to, 
the fourth is the typical match, the fifth the minimax, and 
the final column is the best match. 



m 



O 



On 



typ 



minimax 



Mb 



44.67 0.7987 0.9658 0.9979 0.9974 0.9983 

70.79 0.5458 0.9709 0.9989 0.9989 0.9991 

112.2 0.2308 0.9794 0.9994 0.9993 0.9995 

177.8 0.04864 0.9837 0.9994 0.9993 0.9996 

281.8 -0.06719 0.9893 0.9993 0.9993 0.9996 



for = and a = 0. 

We now analyze the resolution needs of the spinning 
BBH waveforms of which we only have three resolutions 
for each of the spinning configurations in contrast to the 
five available for the non-spinning case. These are labeled 
Acoarso = m/32, Amed = to/38.4, and Ag„o = to/40. 
Later in this paper, we find that it is necessary to in- 
clude the higher modes for the spinning BBH templates; 
and, therefore, as we analyze the spinning BBH wave- 
forms versus resolution we also include its variation with 
0. At inclination angle 6 — only the £ = \m\ = 2 mode 
is present in h, and when ^? 7^ all the modes get mixed. 
We do not include modes £ > 5 in constructing h because 
these were too small to be well resolved for the evolutions 
we use in this paper. The variation of 9 is presented for 
two suggestive cases, 9 = tt/A and 9 ^ tt/2. Values of 9 
are not taken to be exactly on the plane since the radi- 
ation is then linearly polarized and cannot be described 
with two basis vectors and cannot be maximized over the 
phase at that point, but instead at ~ tt/2 = 897r/180. 

For the 9 = case, all the spinning waveforms have a 
-M^minimax > 0.99 ovcr the entire mass range for matches 
between the Acoaisc and Afine and between Amcd and 
Afino. These matches also follow the same trend with 
mass, i.e. decreasing matches with decreasing total mass, 
as reported in Tab. |TT1 For 6* ^ 0, we plot the minimax 
match for the a = 0.8 spin case between both the high- 
est and coarsest and highest and medium resolved wave- 
forms at the two values of 9 in Fig. [51 We only show 
the a = 0.8 waveforms since the matches of the high 
and low resolution of the other spin and non-spinning 
runs are all il/minimax > 0.99 at these angles again over 
the entire mass range. Note, however, that even though 
the matches are very high, they decrease with increasing 
spin. This is expected since we are keeping the number 
of points across the black hole fixed, but as the spin in- 
creases, the horizon area decreases, hence the effective 
numerical resolution of a high initial spin waveform is 
lower than that of a low initial spin waveform. The case of 
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9 ^ -Kjl gave us the lowest matches of any angles between 
and 7r/2 and shows that the minimax match of the 
A/j„e and A^ed resolutions is Mminimax > 0.97. While 
the match of the coarsest resolution run at higher masses 
is Mminimax ~ 0.92, which indicates that the coarse reso- 
lution is not good enough for meeting the threshold when 
a = 0.8 and Q ^ tt/2. For the rest of the paper, we will 
use the finest resolution for all spin cases to ensure that 
the errors due to resolution are no greater 3% in the most 
difficult case to resolve: the high spin, high angle case, 
and should be much lesser in other cases. 
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FIG. 8: We plot the minimax match between the waveforms 
of a = 0.8 for two cases of inclination angle, 9 — n/A and 
9 — n/2 and two different resolutions Afino with Amcd and 
Afine with Acoarse- Corresponding matches for lower spins 
are all > 0.99 and are not shown here. 



Fig. [8] shows two features of the minimax match, 
namely it decreases with increasing 6 and with increas- 
ing TO. The decrease of the match in Fig. [5] at higher 
mass is in contrast to Fig. [7l There are two differences 
between the figures, one is the phase optimization which, 
from Tab. nil we can see does not change the trend of the 
match versus mass. The second difference is the inclusion 
of higher modes. The overall match will be reduced with 
non-zero 9 since the frequencies of the modes increase 
with £, placing stronger resolution requirements. The 
highest spin we analyze, a = 0.8, has the lowest minimax 
matches of all the spin cases, as low as Mminimax = 0.92 
at the largest angle. Large spins demand more resolu- 
tion in general and this is seen particularly when resolv- 
ing modes with £ > 2. For the lower spin runs, a < 0.6, 
the waveforms are well resolved and the matches are all 
> 0.99 over the entire mass range even for 9 ^ tt/2. The 
decrease in match is seen most strongly at the larger mass 
range since this mass range targets the merger and ring- 
down part of the signal. We expect these match values 
and all the matches reported in this section to change 
for waveforms with more gravitational wave cycles since 



the numerical errors will grow as the length of the run 
increases and the accumulating phase errors will likely 
play a larger role. 



C. The Variation of Spinning Templates with 
Inclination Angle 

An observer on the orbital axis of a merging binary, at 
9 = (f) = 0, will see a fully circularly polarized waveform 
where £ = \m\ = 2 is the only nonzero mode. An ob- 
server, however, in the orbital plane of the binary or the 
equatorial plane of the final black hole will see linearly po- 
larized radiation with the contribution from ^ = |to| = 2 
at its minimum. All other observers with intermediate 
orientations would sec elliptically polarized waves that 
are combinations of all the modes. In practice, numeri- 
cal relativity can only resolve a finite set of modes. Since 
the orientation of any given binary will be unknown a. 
priori, we analyze the templates' dependence on the in- 
clination and azimuthal angles. The location of the de- 
tector in the center of mass frame fixes the mode content 
and the orientation of the detector fixes the initial phase 
of the waveform. 

The resolution will now be fixed to the finest for the 
spin cases and the corresponding resolution for the non- 
spinning case. We use the minimax match as defined in 
Eq. pT]) . to see how different the £ ~ \m\ — 2 waveforms, 
given by h{9 = 0, = 0), arc from the full waveforms, 
h{9,(j)), for a given a. Our focus is on 9 since it causes 
a larger variation in the match than (f>. To reduce the 
number of parameters considered, and we set (j) = tt/2 
for the rest of the paper when 9^0. If we were to relax 
this condition, the minimax matches would vary on the 
order of a percent. To observe the variation with 9, we 
first fix the total mass for both templates, hi and to 
IOOA/q. In Fig. [5] we present the minimax matches as a 
fimction of o for templates at 9 = {0, tt/4, tt/S, 7r/2}. By 
holding the target template to 6* = for each spin case, 
we can determine at what angle and spin the matches 
drop below the threshold. In the case oi 9 — 0, the 
templates have the same parameters and Mminimax — 1 
as it should. For a — 0, all the templates match the 9 = 
target within Mminimax = 0.98. The plot also indicates 
that = is close to the full waveform for 9 < tt/3 for 
all of the spin cases. For 9 > 7r/3, however, the match 
drops below 0.98 for a > 0.2. As the inclination with the 
axis increases, the higher modes become important with 
increasing spin. 

To study inclination variation as a function of mass, 
we add the scale with total mass coming from the LIGO 
noise curve. We limit the presentation to two spins, a = 
0.2 and a = 0.8 for clarity. The best, minimax and typical 
matches are plotted as a function of m in Fig. [TO] and 
Fig. [TT] The matches are between hi{9 = 0) and /i2(^i) 
for three inclination angles each 9i = 7r/4, tt/S, ~ 7r/2. 
We plot all three matches to demonstrate the range that 
the matches can take depending on the choice of phase 
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FIG. 9: The minimax match, Mminimax[/ii (9 — 0)], h2{0i)], vs 
a is plotted for the values Oi = {0, 7r/4, n/3, tv/2} for the case 
m = lOOM©. 



optimization. The typical match is given by the lines and 
the minimax and best matches are specified as the lower 
and upper error bars respectively. One can see that the 
minimax match between the two templates sets a lower 
bound on the phase optimized matches. The minimax 
match dips below 0.98 for 6 > tt/3 and m > 100© for 
both spin cases; and, unfortunately, dips below 0.98 for 
the entire mass range for 6* ^ 7r/2 when a — 0.8. 

An even more stringent test is the best match. If the 
best match of a given angle waveform is smaller than 
some threshold then we would know that the higher 
modes are significant and need to be used in creating 
the template bank. We note that for both spin cases the 
best match is below 0.98 for 9 ~ 7r/2 at masses greater 
than I7OM0 indicating the need to include higher modes 
to make a detection at large inclination and high masses. 

One can see from the plots that the lower bound of the 
match between the circularly polarized and the highly 
elliptically polarized waveform is lower in the high mass 
cases. The decrease in the match at larger total masses 
may be indicating that the ringdown is more sensitive 
to the presence of higher modes than the merger itself, 
since the deviation from 0.98 is larger than expected from 
numerical error alone. These figures also suggest that the 
£ — \m\ —2 mode is closer to the full waveform at lower 
masses and lower spins. The fraction of events lost by 
not including the higher modes needs to be calculated to 
predict the full impact. 

The waveforms that we are using have been extracted 
from the numerical solution of the Einstein equation at a 
finite radius. In order to assess how the extraction radius, 
Text, affects the waveforms at a large extraction radius, 
we tabulate the Mminimax of the a = 0.8 waveforms in 
Tab. mil The matches are computed for a template fixed 
at 6* = and the target template at four values of 9 to 



FIG. 10: The range of phase optimized matches between the 
9 = waveform and the full waveform for a low spin case, a = 
0.2. The curves denote typical match Mtyp [/ii (0 = 0),/i2+(^i)] 
for 9i = {7r/4, 7r/3, ~ 7i"/2}. The lower end of the error bar is 
given by the minimax match and the higher end is given by 
the best match. The phase optimization here is done over the 
phase of the hi template. 

measure how much of the variation with 9 depends on 
Text- The variation with the extraction radius is of the 
order of ~ 1.3% in the worst case. There is an interplay 
between better dynamical quality and increased numeri- 
cal error as Text is increased, which is more noticeable in 
modes higher than £ — 2. For the rest of the paper we 
use the waveforms at Text = 30m. 

TABLE III: Minimax matches for m = lOOi\f0 and a = 0.8 
between two waveforms at the same extraction radius but dif- 
ferent inclination angle given by Mminimax[/ii = 0), h2{92)] . 



Text 


6*2 = 


02 = 7r/4 


02 = Tv/3 


e»2 ~ 7r/2 


30 


0.9779 


0.9602 


0.9538 


0.9518 


40 


0.9773 


0.9530 


0.9465 


0.9519 


50 


0.9756 


0.9473 


0.9423 


0.9546 


60 


0.9734 


0.9445 


0.9409 


0.9564 



D. Comparing the Various Spin Configuration 
Waveforms 

Given the computational cost of producing templates 
for BBHs, we investigate how well a reduced template 
matches a target template of arbitrary spin and inclina- 
tion angle at a fixed resolution (A = rn/40). We will 
study several cases, including spin configurations rela- 
tive to each other at different inclination angles for a 
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FIG. 11: The range of phase optimized matches between the 
^ = waveform and the fuU waveform for a high spin case, a = 
0.8. The curves denote typical match Mtyp [hi (6 = 0), /12+ {9i)] 
for 9i — {vr/4, 7r/3, ~ 7r/2}. The lower end of the error bar is 
given by the minimax match and the higher end is given by 
the best match. The phase optimization here is done over the 
phase of the hi template. 



fixed mass of 1OOM0 and then we specialize to a few an- 
gles and study the typical, minimax and the best phase 
optimized matches. 

The simplest approach to compare the spinning wave- 
forms is to fix 6* = for both the target and template. We 
then choose the template to have some spin and calculate 
the match of that template with a target that varies with 
a, such that we calculate MaiinimiiK[hi{ai),h2{aj)] where 
Oi and aj run over all the combinations of spin. The 
minimax matches between all combinations of the spins 
are > 0.995 over the mass range considered. This reflects 
the fact that the binaries approach the same final black 
hole and that the initial spin gets radiated away in modes 
other than £ = \m\ = 2, where they are almost identical 
during the late stages of the binary merger (sometimes 
called universality 42]). This would mean that we could 
not distinguish between the non-spinning and spinning 
waveforms making the non-spinning case sufficient for a 
detection template bank but potentially making param- 
eter estimation problematic. This degeneracy with spin, 
however, does not hold when we include more radiation 
modes. 

We explore the matches between different spin tem- 
plates at different inclination angles, considering first the 
templates with a fixed mass of IOOMq. The variation of 
the minimax match between a template of a = and tar- 
get templates that vary with a is presented in Fig. [T^ and 
with both the template and target varying with a 7^ in 
Fig. [T21 In terms of the match. Fig. [T^ corresponds to a 
minimax match, Afminimax[^i(a = 0), /i2(ai)], and Fig.[T51 
between six combinations of the four spinning configura- 



tions, Mminimax[^i(ai), ^2(aj)] both versus 9. In both 
figures, at 6* < 7r/3 the template is within a Afminimax 
of 0.98 with the targets, making them indistinguishable. 
When the template is fixed at a = in Fig.[T^ the match 
is > 0.98 for all 6 except with a > 0.6, which shows 
a stronger drop-off with angle for higher initial spins. 
Fig. [13] shows that the matches between adjacent tem- 
plates among the a = {0.2,0.4,0.6,0.8} group are better 
than 0.99, while the matches that have a spin param- 
eter difference ai — 02 > 0.4 drop below 0.98 at higher 
angles. These templates are based on short waveforms of 
just two and a half cycles and the matches will likely get 
worse when more cycles are included. We conclude from 
these plots that the £ — \m\ —2 mode continues to dom- 
inate the other modes at small angles, 9 < tt/A, leading 
to similar matches between a spinning and non-spinning 
configuration. As the inclination angle grows, however, 
the presence of higher modes becomes more pronounced, 
breaking the degeneracy between the waveforms from dif- 
ferent initial spin configurations. 
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FIG. 12: The minimax match Mniinimax[^i (a — 0),/i2(ai)] 
as a function of 9 between the waveforms of different spins at 
the same resolution m/40, with the non-spinning waveform at 
resolution m/38.4. The mass of the final black hole for all the 
cases is IOOA/q and the spins run over ai = {0.2, 0.4, 0.6, 0.8}. 
Note the monotonic decrease of the match with angle. 

For completeness, we present the table of the range 
of values over which the phase optimized matches of the 
different spinning waveforms can vary in Tab. IIVI These 
matches are evaluated at 61 7r/2 to report the widest 
range of matches calculated. The table indicates that the 
range of matches is higher at lower masses and tends to 
decrease with mass. 

To place the minimax matches between the spins in 
context with the typical and the best matches, we present 
the typical match Mtyp[/i+i(a = 0), ft,2(a = Oj)] for aj = 
{0.2, 0.8} and 9 = {7r/2, 7r/4} in Fig. [ll The error bars 
are given by the minimax and the best match in order 
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FIG. 13: The minimax match between the waveforms of dif- 
ferent nonzero initial spins, Mininimax[fti(ai)i '*2(aj)] as a func- 
tion of 6. The mass of the final black hole for all the cases 
is IOOM0 and {ai,aj} runs over {0.2,0.4,0.6,0.8}. Note the 
monotonic decrease of the minimax match with the angle 9. 
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FIG. 14: We plot the range and the typical values of phase 
optimized match between the waveforms of different initial 
spins, Mtyp[/i+i(a = Q,9j),h2{ai,6-j)\ as a function of mass 
where Oi = {0.2,0.8} and Oj = {7r/4,7r/2}. 



TABLE IV: The range of phase optimized matches reported 
as (Mininimax, Mbest) between two templates of various spin 
pairs (ai, 02) where both templates have 9 ^ -n 12. 



(m, a2) 



50Mcr 



120AfQ 



300A/( 







(0.0,0.2) 
(0.0,0.4) 
(0.0,0.6) 
(0.0,0.8) 
(0.2,0.4) 
(0.2,0.6) 
(0.2,0.8) 
(0.4,0.6) 
(0.4,0.8) 
(0.6,0.8) 



(0.8940,0.9513) 
(0.8570, 0.9507) 
(0.8595, 0.9530) 
(0.7325, 0.9196) 
(0.9860, 0.9964) 
(0.9742, 0.9927) 
(0.9526, 0.9832) 
(0.9973, 0.9993) 
(0.9808, 0.9956) 
(0.9911,0.9982) 



(0.9228,0.9618) 
(0.8882, 0.9604) 
(0.8898, 0.9618) 
(0.7896, 0.9326) 
(0.9884, 0.9967) 
(0.9751,0.9927,) 
(0.9497, 0.9793, ) 
(0.9959, 0.9988, ) 
(0.9732, 0.9902, ) 
(0.9880, 0.9950, ) 



(0.9446, 0.9701) 
(0.9141,0.9684) 
(0.8234, 0.9427) 
(0.8369, 0.9449) 
(0.9858, 0.9954) 
(0.9665,0.9912) 
(0.9194,0.9681) 
(0.9945,0.9971) 
(0.9611,0.9738) 
(0.9824, 0.9867) 



to demonstrate their variation with mass. As seen in 
Fig. [HI the difference between the minimax and best 
matches is slight at small values of B and a, while for 
large angles it can vary over the order of ^ 3%. The 
variation is much higher at lower masses because there 
are more cycles in the sensitive part of the LIGO noise 
curve making the match more sensitive to the relative 
phasing of the templates. The figure indicates that using 
a non-spinning waveform as the template in detection 
would cause the largest losses when the target has large 
spin and angle. 

We have seen the variation of the match between tem- 



plates of the same spin parameters and different angles 
and different spin parameters and the same angle. In 
setting up a template bank using these numerical wave- 
forms, it may be useful to know if just a couple of tem- 
plates at some specific inclination and spins can cover the 
whole set of waveforms considered here. 

For this study we calculate the minimax matches be- 
tween a template of a given value of spin and inclination 
angle and a target with a spin parameter varying over 
a, = {0.0, 0.2, 0.4, 0.6, 0.8} and a fixed, but different incli- 
nation angle. In Fig. [15] we plot the variation of the min- 
imax match for a template with a\ = 0.0 and Q\ = 7r/4 
with target templates of varying spin and 02 7r/2 as a 
fimction of total mass. We made the choice of the tem- 
plate to have d\ — ii/A since we have concluded that a 
choice of = is going to fail to match well over a range 
of target inclination angles, especially at lower masses. 
We choose the target template to have ^2 ~ 7''/2 to cal- 
culate the worse case scenario for the inclination angle 
of the target. In Fig. [TH] we substitute a template of 
fli = 0.4 and keep all the other parameters the same as 
in Fig. [151 This improves the matches at smaller total 
mass. One can see from the figures that the higher the 
match threshold, the smaller the mass range covered by 
the template. A coarse template bank would do better 
with a — 0.4 and 6 = ■k/A than a = at the same 6 = n/A 
especially for the lower masses and hints that there may 
be an optimal way to lay out this template space. An- 
other feature is that for higher masses, the matches for 
both the cases look similar, as the late-merger/ringdown 
has some universality-independence from the initial spin. 

The use of these axial-spin waveforms as templates for 
detection has some limitations. The waveforms investi- 
gated had the angle between the initial spins kept fixed. 
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FIG. 15: We plot the minimax match between the target 
templates of different initial spins at a fixed 9 with a template 
of a = 0.0, Mmi„imax[/ii(a = 0,6* = 7r/4), /i2(ai, = n/2)] as 
a function of mass. Notice the low matches at the smaller 
masses. 




Q gg I I I I I I I I I I I I I I I I I I I I I I I I 

50 100 150 200 250 

m/M„ 



FIG. 16: We plot the minimax match between the target 
templates of different initial spins at a fixed d with a template 
of a = 0.4, Afminimax[/ii(a = 0.4, = 7r/4), h2{ai, = 7r/2)] as 
a function of mass. 



and it is not clear how the conclusions we draw from this 
study hold as this condition is relaxed. Although it is 
likely that the higher radiation modes will play an even 
bigger role than they do for the relatively simple case 
investigated here. The major limitation however, is the 
small number of gravitational wave cycles present in the 
waveforms. The ISCO frequency is about 0.02/m. The 
mass at which the ISCO frequency hits the seismic wall 
of LIGO noise curve of 40 Hz is about IOOM0. If one 
allows the lower limit in the match integrals to vary or 



be fixed at 0.02/m = 70 Hz for example, then wc can 
trust these matches up to « 60 Mq assuming that we 
will search for events that arc merging in the LIGO band 
with little to no inspiral. Although the results will not 
change too much at high masses even if inspiral is in- 
cluded in the numerical templates, low mass results are 
expected to have lower matches. Finally, another limi- 
tation of these matches is that it assumes identical total 
mass for the template and the waveform being compared 
to it. This assumption needs to be relaxed to estimate 
the "effectualness" of templates for detection, which will 
be done in a future work. 



IV. CONCLUSION 

The parameter space of the BBH system is rapidly 
being explored by the numerical relativity community. 
Given how much time and effort in takes to evolve the 
BBH spacetimes and then to use those waveforms in de- 
tection schemes, it is becoming increasingly important 
to probe the parameter space of the BBH system as the 
waveforms become available. With this motivation, we 
have investigated a series of waveforms from the final 
two to three orbits of a binary-black-hole system made 
up of equal but spinning black holes with matched filter- 
ing. Although the waveforms we use are not long enough 
in terms of gravitational wave cycles to cover the entire 
LIGO frequency band, they shed light on what weight 
different parameters will have in the matched filtering 
procedure using LIGO's initial noise curve. We study 
the dependence of the match on the intrinsic parameters 
of numerical resolution, inclination and spin as well as 
the extrinsic parameters of arrival time and phase. 

Given the numerical convergence properties of our 
waveforms, wc first predict and show that the match has 
a dependence on the resolution of A'^'p, in which a fourth 
order evolution scheme has p = 4. For the resolutions of 
our BBH evolutions, we get very high matches (> 0.98) 
for resolutions above m/32 when a = and above m/38 
for waveforms with a = 0.8. The resolutions are only 
appropriate for the cases we have investigated. For the 
waveforms generated out of the £ = \m.\ = 2 mode, we 
find that the match values between the same waveforms 
at different resolution decreases as the total mass de- 
creases. We only report matches from 50Mq — 3OOM0, 
where 50Mq is set by our initial orbital velocity. 

We then include the variation of the match with the 
inclination angle. It is often convenient and instructive 
to extract the gravitational waveform, '^4, decomposed 
into spin-weighted spherical harmonics as a function of £ 
and TO. Including modes higher than the dominant mode, 
^ = |to| = 2, is equivalent to allowing the inclination an- 
gle of the binary system to vary. Only at 9 = would the 
£ = \m\ = 2 mode be the only harmonic in the waveform. 
The target waveform will have an unknown orientation; 
and, therefore we vary the inclination to study what effect 
this may have on detection. As we include all the modes 
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in the analysis the match decreases, to as low as w 0.85, 
must notably for larger masses and higher spins. Some 
decrease in match is expected since the higher modes re- 
quire more resolution; however, the decrease was greater 
than that accounted for by the resolution alone. 

Finally, in light of the computational effort involved in 
searching over large, densely populated template banks, 
we calculate the match between waveforms of different 
initial spin to see if a reduced set of spinning wave- 
forms would be good enough. The waveforms we use 
in this analysis evolve toward the same Kerr black hole 
(a « 0.66) even though the initial spins of the individual 
black holes vary from 0.0 to 0.8 because of the fixed ini- 
tial ADM angular momentum. This results in very high 
matches, Mminimax > 0.99, between all the spinning and 
non-spinning waveforms for the dominant mode. With- 
out adding information from the higher radiation modes, 
all the spinning BBH waveforms appeared the same in 
terms of matched filtering. 

Once h was constructed using modes £ < 5, this de- 
generacy broke. We constructed two "typical" cases, one 
in which the template had a fixed 6 = -k/2 and a = and 
the other a = 0.4. The templates were matched against a 
target template also at 6 = tt/2 for each spinning wave- 
forms. When a = 0, the matches were as low as 0.88 
versus m, but were > 0.95 for m > IOOMq. By choos- 
ing an intermediate spin of a = 0.4, the matches were 
improved for the lower masses, m < IOOMq, increasing 
the lowest value from 0.88 to 0.96. We speculate that 
this indicates that there may be an optimal layout for 
the template space for a coarse search over spins. The 



matches were always at their highest at IOOMq, although 
not necessarily greater than 0.98. 

In summary, we have investigated the sensitivity of 
matched filtering to the presence of spin and radiation 
modes in a limited case of spinning BBH waveforms. We 
find that if we ignore the modes higher than £ = 2, there 
is very little difference between the spinning waveforms, 
especially at large m due to our choice of initial data. We 
do find, however, that the inclusion of larger i changes 
this picture in two ways. Two templates with the same 
spin but different inclination only match for low values of 
the inclination angle and two templates of the same incli- 
nation angle but different spin only match when the spins 
are within 0.2 of each other. Matching well means that 
Mminimax > 0.98. In general, the higher modes reduce 
the matches more at large masses than at small masses. 
It remains to be seen how these results change when the 
spins are no longer fixed to axial configurations, but more 
elaborate initial data will likely vary more with £. 
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